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ABSTRACT 

We examine the clustering of galaxies around a sample of 20 luminous low redshift (z < 0.30) quasars 
observed with the Wide Field Camera-2 on the Hubble Space Telescope ( HST ). The HST resolution 
makes possible galaxy identification brighter than V = 24.5 and as close as 1" or 2" to the quasar. We find 
a significant enhancement of galaxies within a projected separation of < 100 h^ 1 kpc of the quasars. If we 
model the QSO/galaxy correlation function as a power law with a slope given by the galaxy/galaxy 
correlation function, we find that the ratio of the QSO/galaxy to galaxy/galaxy correlation functions is 
3.8 ± 0.8. The galaxy counts within r < 15 h~ { kpc of the quasars are too high for the density profile to 
have an appreciable core radius (>100 h~ l kpc). Our results reinforce the idea that low redshift quasars 
are located preferentially in groups of 10- 20 galaxies rather than in rich clusters. We see no significant 
difference in the clustering amplitudes derived from radio-loud and radio-quiet subsamples. 

Subject headings: galaxies: clusters: general quasars: general — radio continuum: galaxies 


1. INTRODUCTION 

Over the last two decades, it has been well established 
that quasars are associated with enhancements in the 
galaxy distribution. Historically, this provided the first 
direct observational evidence that quasars were indeed 
cosmological in origin (Bahcall, Schmidt, & Gunn 1969; 
Bahcall & Bahcall 1970; Gunn 1971 ; Stockton 1978). Over 
the years, considerable evidence has accumulated that low- 
redshift (z < 0.4) quasars reside in small to moderate groups 
of galaxies rather than in rich clusters (cf. Hartwick & 
Schade 1990; Bahcall & Chokshi 1991, and references 
therein). At higher redshifts, there is a marked difference in 
the environments of radio-loud and radio-quiet quasars 
(Yee & Green 1987). At redshifts z > 0.6, radio-loud 
quasars are often found in association with rich clusters 
(Abell richness R > 1), while radio-quiet quasars appear to 
remain in smaller groups, or, perhaps, in the outer regions 
of clusters (Boyle, Shanks, & Yee 1988; Yee 1990). 

The galaxy environment around quasars provides many 
important clues as to what triggers and sustains their 
central engines. As first suggested by Toomre & Toomre 
(1972), mergers and interactions of galaxies can provide an 
efficient mechanism for transporting gas into the inner 
regions of a galaxy or quasar. There have been attempts to 
model the interaction/merger rates of ordinary galaxies in 
order to explain the luminosity function of quasars (De 
Robertis 1985; Roos 1981b; Carlberg 1990), and the rapid 
evolution of the merger/interaction rate in clusters with red- 
shift may provide a natural explanation of the strong evolu- 
tion of clustering observed around radio-loud quasars 
(Stocke & Perrenod 1981; Roos 1981a). Knowledge of the 
galaxy environment around quasars is also important for 
understanding the large-scale distribution of quasars and 


1 Based on observations with the NASA/ESA Hubble Space Telescope , 
obtained at the Space Telescope Science Institute, which is operated by the 
Association of Universities for Research in Astronomy, Inc., under NASA 
contract NAS 5-26555. 


how it relates to the structure seen in galaxy su \ 
(Bahcall & Chokshi 1991). 

In this paper, we examine the galaxy environment ar >i 
20 nearby (z < 0.3) bright quasars observed with the A 
Field and Planetary Camera-2 (WFPC2) of the H d 
Space Telescope (HST). These fields were imaged as p r 
an ongoing investigation into the nature of the 1 
environment of quasars (Bahcall, Kirhakos, & Schr ; ei 
1994, 1995a, 1995b, 1996a). The exceptional resolution 
the HST images allows companion galaxies to be det x 
at very close projected separations, in some cases r ~ / 

kpc (<2"), and galaxy/star separation to be perfo n 
down to V ~ 24.5. The goal of the work presented heri i: 
quantify the excess of galaxies associated with h 
quasars. The outline of this paper is as follows. A b 
description of the quasar sample used in our analy a: 
given in § 2. In § 3.1, we argue that galaxy counts are ii o 
sistent with being drawn from a uniform background 
§ 3.2, we strengthen this conclusion by correcting the ci-u 
for the background contamination. We also presen 
excess galaxy counts above the background in annuli 
projected separation. From these counts, we quantif 
amplitude of the galaxy clustering around the quasar 
§ 3.3 in terms of a quasar-galaxy spatial cross-correl it 
amplitude. We discuss our results and their relation tc i 
vious work in § 4. 

2. DATA 

The sample of objects analyzed in this paper consi: o 
20 of the intrinsically most luminous (M v < —22.9 
H 0 = 100 km s _1 Mpc~\ Q 0 = 1) nearby (z < 0.30) v .c 
quiet and radio-loud quasars selected from the Veron-i \ 
& Veron (1991) catalog. Table 1 lists the individual qu is 
and their redshifts. The quasars span the redshift i n 
0.086 < z < 0.29, with a median redshift of z med = 0 . 
Details of the observations and the resulting images for 
of the 20 fields have been presented in Bahcall et al. ( 9 
1995a, 1995b); the remaining observations will be presen 
in Bahcall, Kirhakos, & Schneider (1996b). Six ol 
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TABLE 1 
Quasar Sample 


Name 

Date 

Exp. (s) 

My 

z 

PG 00523-251 

05 Dec 94 

1400 

-23.0 

0.155 

PHL 909 

17 Oct 94 

1400 

-22.9 

0.171 

NAB 0205 + 02 

26 Oct 94 

1400 

-23.0 

0.155 

0316-346 

20 Nov 94 

1400 

-24.5 

0.265 

PG 0923 + 201 

23 Mar 95 

1400 

-23.0 

0.190 

PG 0953+414 

03 Feb 94 

1100 

-24.1 

0.239 

PKS 1004+13“ 

26 Feb 95 

1400 

-24.3 

0.240 

PG 1012 + 00 

25 Feb 95 

1400 

-23.0 

0.185 

HE 1029-1401 

06 Feb 95 

1400 

-23.2 

0.086 

PG 1116 + 215 

08 Feb 94 

1100 

23.7 

0.177 

PG 1202 + 281 

08 Feb 94 

1100 

-23.0 

0.165 

3C 273* 

05 Jun 94 

1100 

-25.7 

0.158 

PKS 1302-102“ 

09 Jun 94 

1100 

-24.6 

0.286 

PG 1307 + 085 

05 Apr 94 

1400 

-23.1 

0.155 

3C 323.1“ 

09 Jun 94 

1100 

-22.9 

0.266 

PG 1309 + 355 

26 Mar 95 

1400 

-23.2 

0.184 

PG 1402 + 261 

07 Mar 95 

1400 

-23.0 

0.164 

PG 1444 + 407 

27 Jun 94 

1100 

-24.0 

0.267 

PK.S 2135-147* 

15 Aug 94 

1400 

-23.5 

0.200 

PKS 2349-014“ 

18 Sep 94 

1400 

-23.3 

0.173 


“ Radio loud. 


quasars (see Table 1) are radio-loud, while the remainder 
are radio-quiet (Kellermann 1989). 

Each quasar field was imaged with the WFPC2 through 
the F606W filter, which is similar but slightly redder than 
the V bandpass (5940 A, A = 1500 A). The quasars were 
positioned within 4" + 1'.'2 of the center of Wide Field 
Camera CCD 3 (WF3). Simultaneous images were obtained 
in the adjacent CCD chips 2 and 4 (WF2 and WF4, 
respectively), which together with WF3 form an “ L ” 
shaped image (see Fig. 1). Each chip has 800 x 800 pixels 
and an image scale of 0"0996 pixel -1 at the chip’s center; 
this corresponds to spatial resolution of 2.1 h 1 kpc 
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Fig. 1. — Mosaic of the galaxy counts in WF2, WF3, and WF4 ( upper 
lefty lower lefty and right panels, respectively) for all twenty fields. In a given 
field i , we count only those galaxies in the magnitude range m^zj - 1 to 
mjzj + 2, where m^zj is the apparent magnitude of an L + galaxy at the 
redshift of the quasar, z t . 


arcsec - 1 at a redshift z = 0.20 (fi 0 - 1). The effective areas 
(areas not covered by pyramid shadows) of WF2, WF3, and 
WF4 are 1.59, 1.60, and 1.59 arcmin 2 , respectively. More 
detailed information on the WFPC2 and its photometric 
system can be found in Burrows (1994), Trauger et al. 
(1994), and Holtzman et al. (1995a, 1995b). The relatively 
long exposures (1 100 or 1400 s), combined with the excellent 
spatial resolution, allowed galaxies to be identified in the 
images down to limiting magnitude m(F606W) < 24.5 and 
as close as <1" or 2" from the central quasar. We per- 
formed aperture photometry on the field galaxies; circular 
apertures with radii of 0'3-10" were used, as appropriate. 

3. GALAXY COUNTS AROUND QUASARS 
3.1. Raw Counts : Evidence for a Strong Enhancement 
If galaxies are distributed around low-redshift quasars 
with a power-law distribution, £(r) ~ (r/10 h 1 Mpc) 1 77 
(as suggested by, e.g., Yee & Green 1987), then there will be 
a very strong enhancement of the counts within 30" (r ~ 60 
h~ l kpc at 2 = 0.2) of the quasar. Moreover, because the 
centers of WF2 and WF4 are offset from the quasar, there 
will be an enhancement of galaxies in WF3 relative to WF2 
and WF4. The background galaxy counts rise steeply with 
magnitude, and this will dilute any signal of excess galaxies. 
Much of the background contamination can be removed 
simply by counting only those galaxies with apparent mag- 
nitudes in the range that is likely to be physically associated 
with the quasars. A good compromise between eliminating 
too many associated galaxies and minimizing the effect of 
foreground/background interlopers, is to count galaxies in 
each field, i, that are in the magnitude range m^z*) - 1 to 
m *( z i) + 2, where m, )c (z i ) is the apparent magnitude of an L* 
galaxy at the redshift of the quasar, z x . The mean (median) 
value of m*(z,) for our sample is 18.2 (18.1); the total range 
of m + (z f ) is 16.4 < m + (z f ) < 19.2. In computing m # (z), we 
have taken an L * galaxy to correspond to an absolute mag- 
nitude M*( F606W) = -20.75 and have used the K correc- 
tions between absolute and apparent magnitudes given in 
Fukugita, Shimasaku, & Ichikawa (1995). 

Figure 1 shows the positions of all galaxies in our twenty 
fields, with m + (z*) - 1 < m < m^z,*) + 2; geometric correc- 
tions were applied according to Holtzman et al. (1995a). In 
the panel containing the quasar (WF3, lower left), there is a 
clear excess of galaxies (50 galaxies) relative to WF2 (upper 
left) and WF4 (18 and 17 galaxies, respectively). The 
hypothesis that the counts in all three chips are drawn from 
an underlying Poisson distribution with a common mean 
leads to a maximum likelihood mean per chip of 28.4 and 
X 2 = 24.5; the probability that yf for one degree of freedom 
is this large by chance, is only P = 7 x 10 -7 . Thus, without 
any detailed modeling, we can rule out the possibility that 
the counts are random fluctuations in the background dis- 
tribution. In the following sections, we attempt to make this 
conclusion progressively more quantitative, by first model- 
ing the background galaxy distribution, and then intro- 
ducing a model for the spatial distribution of galaxies 
around the quasars. 

3.2. Correcting for the Background Galaxies 

In order to further quantify the excess of galaxy counts 
around the quasars, we need an estimate of the contribution 
from background galaxies. In Figure 2, we show the galaxy 
counts (d 2 N/dmdQ arcsec -2 ) versus magnitude for the 
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Fig. 2. — Number counts, d 2 Nfdmdi 2, versus m(F606W). Filled circles 
represent the counts determined directly from the off chips in each QSO 
field. Poisson errors are shown. The line is a least squares fit to the faint 
counts, with m(F606W) > 21.5. 


“ off quasar” chips WF2 and WF4, in bins Am(F606W) = 
0.5. The counts at w(F606W) > 21.5 mag are well approxi- 
mated by a power law, log 10 d 2 N/dmdQ = — 10.8 + 
0.33m(F606W); the counts at brighter magnitudes are in 
excess of those obtained by extrapolating the faint power- 
law fit. Because the separation of the centers of chips 2 and 4 
is not large (101" corresponds to 213 h~ l kpc at z = 0.2), the 
counts in these chips are a combination of both background 
galaxies and the (relatively bright) galaxies physically 
associated with the quasars, and, hence, they yield an over- 
estimate of the true background. We have compared the 
power-law fit with the counts derived from the HST 
Medium Deep Survey (MDS) (Griffiths et al. 1994; S. Caser- 
tano 1995, private communication), which covers a much 
larger area of the sky. The agreement is good, and in the 
remainder of this paper, we will adopt the power law in 
Figure 2 as our estimate of d 2 N/dm dCl. The agreement with 
the MDS is also a useful consistency check for systematic 
errors in our derived magnitudes. 

In the 20 fields, there are 1 1 galaxies in the range m*(z,) 
— 1 < m < + 2 within a projected separation less 

than 25 h~ l kpc of the quasar; the total number expected 
from the d 2 N/dmdQ power-law relation is only 0.99. The 
probability of the observed counts being a Poisson realiza- 
tion of the background is extremely small, P — 8 x 10“ 9 , 2 
order of magnitudes smaller than our previous estimate, 
obtained by neglecting the background. This is also a much 
stronger constraint than the upper limit given in Bahcall et 
al. (1995a), P = 2 x 10~ 2 , based on eight fields and an 
(over)estimate of the background obtained from the counts 
in WF2 and WF4. 

We have also computed counts within projected separa- 
tions of 5 and 10 /i _1 kpc. In the 20 fields, we find 2/0.039 
and 5/0.16 (number/number expected in the background), 
with r < 5 h~ l kpc and r < 10 h~ 1 kpc, respectively; the 
corresponding random probabilities are P = 7 x 10 4 and 


P = 1 x 10“ 7 . There is a suggestive difference in the c i 
derived from the radio-loud and radio-quiet subsan p 
The radio-quiet quasars (14 fields) had 0/0.027 (P = 0 
and 2/0.11 (P = 5 x 10" 3 ) for the 5 and 10 h 1 kpc cc n 
the corresponding numbers for the radio-loud sa n 
(6 fields) were 2/0.012 (P = 8xl0 5 ) and 3 0 
(P = 2 x 10“ 5 ). 

In Figure 3, we show the fractional excess of gal i. 
above the background, <<5N/jV(r)> = <JV,-(r)/iV^r)) - 
obtained by averaging the counts in the 20 fields in bi i 
15 /i _1 kpc projected separation. Here N t {r) is the a :l 
count of galaxies with projected separations bet v 
r — Ar/2 and r + Ar/2 in the zth field, and N h Jr) i 
expected background contribution. Again, we have o 
considered those galaxies in the apparent magnitude i u 
m*(z,) — 1 to m + (z ( ) -F 2. 

From Figure 3, we can immediately draw two co it 
sions: first, there is a significant excess of galaxies u;t; 
projected separations of r < h ~ i kpc from the quasars i 
second, there appears to be no significant difference i. i 
galaxy counts for the radio-loud and radio-quiet s; 
samples for 10 h~ x kpc < r < 100 h 1 kpc. In the n 
section, we quantify the observed clustering in terms c f ; 
spatial quasar/galaxy cross-correlation function. 


3.3. Estimating the Spatial Clustering Amplitude 

A detailed derivation of the relation between an- u 
counts and a spatial distribution of galaxies in terms o 
cross-correlation function is given in Longair & Sei ii 
(1979). Briefly, one assumes that the galaxy distribi ti 
(above the background) around the quasar is specified b 
quasar/galaxy cross-correlation function, £ qg (r, z). 1 
observed excess in projected separation in the /th fie d 
then obtained by integrating f w (r, z) over the redshift d 
tribution of galaxies (dA r /dz) i in the apparent magni i 
range 0*(z,) - 1 < m < >n*(zj + 2], 


NM/N b Jr) - 1 


\dOdz{d. Vjdz) l ^(s, z) 
j d6 dz(dA r fdz) i 


In this equation, s is the comoving separation betweei 
line of sight, x(z), and the quasar located at coordinate . ( 
s 2 = {(x(z,.) - x(z)) 2 /G(x(z))} + x(z) 2 9 2 ; 0 is the < n 

between the line of sight and the quasar; and G(x) n 
function that describes the degree of spatial curva u 
G(x) = 1 — (H 0 x/c) 2 (Q 0 — 1) (cf. Peebles 1980, eq. 5(1 
The range of integration in 6 corresponds to the angles il 
span the bin of projected separation, i.e., from (r — A 7 
d a (Zi) to (r -h Ar/2)/d a (z f ), where d a (z f ) is the angular dian e 
distance to the quasar. 

The redshift distribution (d~V jdz\ (assuming no gal; \ 
are created or destroyed) is given by 



where </>(z ( ) is the integral of the luminosity function n 
the absolute magnitudes [M *(2 f ) - 1 < M < MJzj - 
corresponding the apparent magnitude range [mj: ) 

1 < m < m^Zj) + 2 ], 

The cosmological model (Q 0 , A, H 0 ) enters in the al u 
equations implicitly in the volume element, angular d 



472 


FISHER ET AL. 


Vol. 468 



Fig. 3.— Average galaxy counts < <$N/N(r)> of the full quasar sample and radio and radio-quiet subsets. Solid points are the counts computed, using a 
power-law fit, to the faint counts of chips 2 and 4 for the background (cf. Fig. 2); the open symbols (shifted 0.025 in the log) show the counts derived using the 
background counts from the Medium Deep Survey. Error bars show the scatter between the different fields. Curves are the predicted correlation functions if 
the galaxy/quasar cross correlation function is a power law with index y = 1.77 (so/td)> & n exponential surface density (dashed), or a modified Hubble profile 
(dotted). All calculations assume Q 0 = 1 and A = 0. 


tance, and luminosity distance. We adopt a canonical 
model, with Q 0 = 1, H 0 — 100 km s 1 Mpc -1 , and A = 0. 
The basic results of our study are largely independent of this 
choice. The redshift distribution is computed using a non- 
evolving Schechter (1976) luminosity function with a faint 
end slope of a = 0.97 (Loveday et al. 1992). The K- 
corrections necessary for computing the relation between 
absolute and apparent magnitudes have been taken from 
Fukugita et al. (1995). In the results that follow, the H 0 
dependence is explicitly indicated. The derived clustering 
amplitudes remain within the stated errors, as Q 0 is varied 
from 0.1 to 1.0. 

The cross-correlation function, <^(r, z\ is usually 
assumed to evolve with redshift. A convenient model for 
this evolution is to assume that on small scales the clus- 
tering is constant in physical coordinates, i.e., the number of 
excess galaxies around the quasar, n(z)£ qg (r, z), is a constant 
[ n(z ) is average number density at the epoch z]. Thus, for 
any assumed shape of the correlation function, F(r), £, qg (r, z) 
evolves as 

z| - (Thf frh] - (4) 


The factor of (1 4- z) 3 compensates for the change in the 
mean number density, while the redshift factor in the argu- 
ment of F(r) is simply a matter of convention; one usually 
specifies the shape in terms of physical separations, while in 
(1) we have specified £ qg (r, z) in terms of comoving separa- 
tion. 

We consider three different models for F(r). The first 
model is a power law, F(r) = B qg r~ 1 77 , with a slope equal 
to that of the galaxy/galaxy correlation function (Davis & 
Peebles 1983) and the amplitude taken as a free parameter. 
The second model is an exponential surface density of gal- 
axies, p(r) = p 0 exp ( — r/r c ), which, after deprojection by a 
standard Abel inversion, corresponds to F(r) = 
p 0 /nr c K 0 (r/r c ) (K 0 [x] being the modified Bessel function). 
This model was proposed by Merrifield & Kent (1989) as a 
typical galaxy profile around centrally dominant cluster 
galaxies. We adopt their best estimate of a core radius of 
r c — 100 h 1 kpc and vary the amplitude p 0 . Finally, we 
consider a modified Hubble profile F(r) = A\\ + 
(r/r c ) 2 ] " 3/2 (Binney & Tremaine 1987, eqs. 2-37), with r c = 

1 00 1 kpc and A as a free parameter. These three models 

have different behavior at small r. Both the power-law and 
exponential models diverge as r 0 [although the later 
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TABLE 2 

Clustering Amplitudes 


Sample (B q9 >* 1 - a Range <B qg )KB qg ) h 


All quasars 75 60-93 3.8 ± 0.8 

Radio loud 84 57-117 4.2 ± 1.5 

Radio quiet 72 53- 92 3.6 ±1.0 


a In units of [h 1 Mpc] 1,77 . 
b <£^> = 19.8(/T l Mpc) 1 - 77 


does so only weakly, K 0 (x) ~ — ln(.x)], while the modified 
Hubble model asymptotically approaches a constant. 

We varied the amplitude of the models in order to 
achieve a maximum likelihood fit to the excess galaxy 
counts in bins of projected separation of 15 h~ 1 kpc. The 
limited number of fields prevented us from varying more 
than one parameter for each model. Figure 3 shows the 
resulting best-fit models; the amplitudes of the power-law 
model are given in Table 2. At separations r > 30 h~ l kpc, 
the excess counts are relatively flat and all three models for 
the correlation function fit the data. However, the counts in 
the innermost bin, r < 15 /i -1 kpc, lie significantly above 
the counts at larger separations; the rise in the counts at 
small radii is particularly striking in the five radio-loud 
fields. The power law is the only model for £ qg (r) considered 
that rises steeply enough to account for this excess. 

4. DISCUSSION 

If the quasars were distributed like typical galaxies, then 
the derived value of would be equal to the amplitude 
of the galaxy/galaxy correlation function, < B qg ) ~ 20 
(Davis & Peebles 1983). A higher value of < B qg ) suggests 
that the quasars lie preferentially in regions of above 
average galaxy density. Following Bahcall & Chokshi 
(1991), we can convert our values of B qq into an estimate of 
the typical richness of quasar galaxy environment. Here 
we define richness as the number of L* galaxies associated 
with the quasar; this is given by a simple integral over the 
correlation function, N = 471/?* J £ qg (r)r 2 dr, where the limits 
of integration are r = 0-1.5 h~ 1 Mpc (the traditional Abell 
radius), and n* « 1.5 x 1CT 2 h 3 Mpc~ 3 is the number 
density of L* galaxies. Using the values of < B qg ) in Table 2, 
we find the quasars typically reside in groups of 16 25 gal- 
axies. These numbers should be compared to the typical 
richness of Abell clusters, which have 30-49 and 50-79 
members for richness classes R ~ 0 and R = 1, respectively. 
Moreover, this estimate is most likely an upper limit, since 
there is evidence that the galaxy profile around quasars falls 
off more steeply than r~ 1,77 for r > 0.25 h~ l Mpc 
(Ellingson, Yee, & Green 1991). In order to test the robust- 
ness of our results to a steepening of the galaxy profile, we 
fitted a double power-law model with slope —1.77 for 
r < 0.25 Mpc and —3 for r > 0.25 h~ l Mpc. The best- 
fit amplitudes, B qgy for this model increased by about 10% 
(still well within the quoted 1 cr errors), yet, the number of 
inferred bright galaxies associated with the quasars 
decreased to eight, owing to the steeper profile at large 
separations. 

The derived amplitudes, (B qg ) y in the pure power-law 
case are somewhat larger than the estimates by Yee & 


Green (1987). They examined the clustering of nine l a 
loud and 16 radio-quiet quasars in the redshift a 
0. 15 < z < 0.30. On scales of ^ 20 500 h~ x kpc, they t 
mated <B W > ^ 60 ± 20 and (B qg ) ^42+14 for r u 
loud and radio-quiet quasars, respectively. Hayman (A 
derived galaxy counts around low-redshift (z < 0.3) qtar 
from the Palomar Sky Survey prints. He found thin 
ratio of the quasar/galaxy and galaxy/galaxy angular co 
lation functions was 3.1 ± 0.6; if the quasars and gam 
have similar selection functions, this translates to an i 
mate similar to Yee & Green of (B qg ) ~ 61 ± 12. Freu d 
Gunn (1983) analyzed a sample 25 low-redshift qu is 
(z < 0.35) selected from 1.2 m Palomar Schmidt plates ; 
they concluded that < B qg ) = 25 + 12; they also ana-y 
the data set of Stockton (1978; 27 quasars with red d 
z < 0.45 selected from the red Sky Survey prints) ; 
derived, via the same analysis, B qg = 19 ± 40. Our mea a 
ments are, with the exception of Yee & Green’s value 
radio-quiet quasars, consistent within the quoted esr 
The slightly higher clustering amplitude we derive fc 
radio-quiet subsample may be a result of our sample U 
the subset of the most luminous quasars. 

Yee & Green (1987) found that the clustering ampl i 
of galaxies around radio-loud quasars increased by a f a 
~3 between z ~ 0.4- 0.6, and at z ~ 0.6, radio-loud qu o 
are found in environments as rich as Abell class R - 
Optical quasars do not evolve as rapidly (Boyle et al. 1 >; 
perhaps indicating a different formation scenario. I I 
been suggested that quasars and active galactic nuclei n 
be triggered by interactions (e.g., Toomre & Toomre 9 
Stocke & Perrenod 1981 ; Roos 1981a, Yee 1987). This < 11 
a simple explanation for why the quasars are typicall i 
found in rich clusters at low redshifts; the high-vel u 
dispersion of such clusters leads to a low-interaction ra c 

The HST WFPC2 is an excellent instrument for exie 
ing the present analysis to fainter, low-redshift, quas. 
This extension would improve the counting statistics , i 
also providing information regarding possible correla u 
of the quasar environment with luminosity. The imagi q 
moderate redshift quasars (z < 0.6) could be accompLsl 
by HST with single orbit exposures. This imaging uo 
provide a more direct comparison with previous grcu 
based work, and it would increase our knowledge oi 
evolutionary history of the quasar environment, 
increased knowledge of the quasar environment wou 1 
useful in using quasar observations to probe large- x 
structure at higher redshifts. 
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